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Section  1 

BACKGROUND  AND  INTRODUCTION 

1.1  BACKGROUND 

Many  of  the  initial  studies  concerning  the  response  of  structures  subject  to 
high  intensity  heating  assumed  material  removal  and/or  burnthrough  to  be  the 
primary  structure  damage  mechanism.  In  reality,  structural  failure  due  to 
the  reduction  of  material  allowables  at  elevated  temperature  and  the  resultant 
stress  redistribution  are  more  probable  and  can  result  from  lower  intensity 
heating.  In  more  recent  investigations,  (References  1  and  2)  analytical  and 
experimental  studies  were  conducted  to  evaluate  failure  of  simple  plates  and 
built-up  beams  due  to  high  intensity  heating.  The  primary  objectives  of  these 
studies  were  to  investigate  failure  mechanisms  and  demonstrate  the  ability  of 
computer  codes  to  perform  rapid  heating  analysis. 

The  methodology  was  applied  in  the  vulnerability  assessment  of  various  missiles. 
These  studies  showed  that  catastrophic  structural  failure  could  be  produced 
by  heating  the  structure  to  a  temperature  much  lower  than  that  required  to  cause 
melt  or  burnthrough.  Various  areas  on  missiles  which  are  susceptible  to  such 
damage  include  the  rocket  motor  case,  elevons,  wings  and  guidance  sections. 

A  plan  was  formulated  to  perform  ground  based  tests  on  the  various  structural 
components  to  determine  their  vulnerability.  In  particular,  simulation  tests 
were  performed  on  instrumented  pressurized  and  unpressurized  rocket  motor  cases 
using  a  graphite  heater  radiation  source. 

Analysis  to  predict  the  response  of  the  rocket  motor  cases  and  to  assess  the 
sensitivity  of  the  predicted  response  to  various  heat  sources  was  recently  per¬ 
formed  (Reference  3).  An  interest  in  quantifying  the  data  and  modeling 
requirements  necessary  for  accurate  analysis/experiment  correlation  played  a 
key  role  in  formulating  the  present  study.  The  present  investigation  addresses 
the  role  and  appropriate  form  of  material  property  data  to  be  used  for 
structural  response  and  failure  prediction  due  to  short  time  heating. 
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1.2  INTRODUCTION 

Prediction  of  the  failure  of  structures  which  are  subjected  to  a  combination 
of  1)  mechanical  loads  and  2)  thermal  loads  due  to  rapid  (1-10  seconds),  high 
Intensity  heating  requires  knowledge  of  load  and  temperature  histories  and  a 
thorough  understanding  of  the  response  of  the  material  to  these  load  and 
temperature  histories.  For  the  problem  of  rapid  heating  of  flight  structures 
such  as  missiles  and  aircraft,  structural  failures  typically  occur  at  temper¬ 
atures  between  50-100%  of  the  material  melt  temperature.  These  strength  or 
stiffness  type  failures  are  caused  by  a  load  redistribution  and  reduced  load 
carrying  capability  due  to  the  degradation  of  structural  properties  of  the 
material  at  high  temperature,  an  important  consideration  even  in  the  short 
heating  time. 

Once  the  loads,  geometry  and  materials  are  known,  it  is  up  to  the  structural 
analyst  to  choose  the  appropriate  code  or  analysis  technique  to  predict  the 
response  of  the  structure.  For  failure  prediction, the  analyses,  in  general, 
must  be  capable  of  handling  nonlinear  temperature  and  time  dependent  material 
behavior  and  large  displacement  effects. 

The  ability  of  a  code  to  perform  the  required  analysis  can  be  established  by 
comparing  code  predictions  with  the  results  from  experiments.  Satisfactory 
correlation  depends  upon  essentially  three  factors: 

1.  Accurate  experimental  structural  response  data. 

2.  Suitable  code  or  analysis  capability. 

3.  "Appropriate"  material  property  data  for  code  input. 

Previous  studies  (Reference  1)  have  emphasized  the  predictive  ability  of  the 
structural  response  codes  using  available  material  property  data.  The 
objective  of  this  study  Is  to  1)  Identify  a  possible  need  for  further  high 
temperature  material  property  data  for  short  time  heating,  2)  define  an 
appropriate  format  for  the  data  for  use  in  the  codes,  3)  show  how  this  data 
can  be  obtained,  and  4)  evaluate  the  degree  of  improvement  in  the  response 
prediction  that  can  be  achieved  using  these  new  data. 
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1.2.1  Material  Properties 

It  Is  well  known  that  particular  material  properties  can  greatly  Influence  the 
stress  and  deformation  states  in  a  structure  just  prior  to  failure.  In  indeter¬ 
minant  structures,  these  failures  often  occur  outside  of  the  intensely  heated 
region.  A  large  uncertainty  about  the  appropriate  material  property  data 
still  exists  In  the  mind  of  the  analyst  despite  the  importance  of  this  data. 
There  are  at  least  three  basic  reasons  for  the  lack  of  appropriate  data. 

1.  Most  material  property  data  have  been  obtained  for  use  in  sure  safe 
design  and  analysis.  This  means  that  data  for  structural  materials  are 
available  only  at  the  operating  temperatures  of  standard  structures. 

2.  Data  In  the  literature  are  generally  presented  In  a  form  that  Is 
appropriate  for  linear  analyses  that  are  typically  employed  for  sure  safe 
design.  For  example,  only  the  elastic  moduli  and  yield  stress  are  given  as  a 
function  of  temperature.  For  failure  prediction,  it  is  necessary  to  account 
for  the  full  range  of  stress-strain  behavior  at  temperature. 

3.  Data  in  the  literature  do  not  account  for  time  dependence  in  the  short 
time  region  of  interest.  Most  data  of  this  type  are  for  long  service  times 
appropriate  for  aircraft.  This  problem  is  closely  related  to  the  lack  of 

high  temperature  data  where  short  time  creep  effects  become  important. 

The  current  emphasis  is  placed  on  the  importance  of  mechanical  properties.  It 
should  be  noted  that  certain  physical  properties  such  as  specific  heat, 
conductivity  and  thermal  expansion  are  also  important  in  the  prediction  of 
temperature  within  the  structure  and  thermal  deformation.  In  some  cases,  these 
data  are  also  not  available  at  high  temperature;  however,  the  type  of 
("appropriate")  physical  property  data  necessary  for  the  analysis  is  not  in 
question. 

An  important  task  in  this  study  will  be  to  establish  the  "appropriate"  mechani¬ 
cal  material  property  data.  By  "appropriate"  we  mean  that  data  which  when 
used  in  a  theoretically  valid  constitutive  model  in  existing  codes  will  allow 
one  to  accurately  predict  the  response  of  the  structure.  In  this  context, 
material  property  data  are  assumed  to  be  separable  into  three  types: 

1.  time  independent  -  temperature  and  load  Independent  strain  (elastic/ 
plastic  ) 


2.  time  dependent  -  temperature  and  load  dependent  strain  (creep) 

3.  temperature  dependent  free  expansion  strain. 

The  stress-strain-temperature-time  data  are  generally  determined  from  simple 
one-dimensional  experiments  and  the  appropriate  assumptions  based  on  physical 
Insight  and  mathematical  necessities  are  made  to  extend  these  data  to  multi - 
axial  stress  states.  This  study  will  focus  on  one-dimensional  behavior  both 
experimentally  and  analytically  as  a  springboard  to  understanding  multi¬ 
dimensional  behavior. 

1.2.2  Modeling  and  Computer  Codes 

Several  general  purpose  and  special  purpose  structural  analysis  codes  exist 
which  have  the  theoretical  capability  of  performing  the  various  elements  of 
a  rapid  heating  structural  response  analysis.  Two  general  purpose  codes 
(MARC  and  ANSYS)  were  selected  as  being  representative  of  codes  which  have 
the  capability  of  performing  the  complete  time- temperature-load  analysis  of 
a  complex  structure  given  the  proper  input  data.  These  codes  along  with  a 
special  purpose  one-dlmens lonal  code  (CREEPARHS-creep-elastlc-plasilc  analysis 
of  rapidly  heated  specimens)  developed  at  MDAC  will  be  used  to 
predict  the  response  of  a  one  dimensional  bar  subjected  to  time  varying 
load  and  temperature  profiles  typical  of  those  seen  by  flight  vehicle 
structures. 

1.2.3  Program  Plan  and  Report  Organization 

This  study  will  establish  and  demonstrate  a  systematic  procedure  which  can  be 
used  to  assess  the  suitability  of  mechanical  materiel  property  data  as  used 
by  standard  structural  response  codes  to  predict  the  response  of  complex 
structures  subject  to  rapid  heating.  The  procedure  begins  with  ar.  assessment 
of  the  ability  of  constitutive  models  In  the  code  to  simulate  the  actual 
material  response,  e.g.,  time  dependent  elastic-plastic-creep  In  stress/ 
temperature  regimes  where  creep  effects  arc  Important.  Once  the  "appropriate" 
material  property  data  requirements  are  defined,  specimens  which  will  provide 
the  required  data  are  designed  and  the  "appropriate"  tests  performed.  The 
ability  of  this  data/code  combination  to  accurately  predict  structural 
response  Is  then  demonstrated  by  comparing  the  analytical  results  with 
results  obtained  from  simple  simulation  experiments  for  typical  structural 
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load,  temperature  and  time  histories.  The  analysis  and  experimental  correla¬ 
tion  on  the  one-dimensional  specimens  serve  as  a  basic  check  on  the  material 
models,  code  operation  and  experimental  results.  One  can  then  precede  with 
some  confidence  to  analyze  more  complex  structures.  A  flow  diagram  of  study 
tasks  Is  shown  In  Figure  1.  The  following  sections  will  present  the  results 
of  the  study.  A  summary  of  the  work  Including  conclusions  and  recontnendatlons 
Is  given  In  Section  2.  A  survey  of  available  analytical  methods  and  models  Is 
presented  In  Section  3  followed  by  a  description  of  the  test  equipment  and 
data  reduction  scheme  In  Section  4.  The  material  property  and  simulation 
experiments  are  described  In  Sections  5  and  6  and  a  discussion  of  the 
correlation  between  analysis  and  experiment  Is  given  In  Section  7. 


Figure  1.  Flow  Diagram  of  Program  Tasks. 
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Section  2 

SUMMARY,  CONCLUSIONS,  RECOMMENDATIONS 


2.1  SUMMARY 

An  Investigation  of  the  theoretical  requirements  for  an  accurate  response 
analysis  of  a  structure  which  Is  subject  to  a  combination  of  mechanical 
and  thermal  loads  due  to  rapid  heating  was  undertaken.  An  assessment  of 
the  ability  of  current  nonlinear  structural  codes  (MARC,  ANSYS)  to  predict  the 
aforementioned  response  was  made  In  light  of  the  theoretical  requirements.  A 
critical  evaluation  of  the  suitability  of  existing  mechanical  material  pro¬ 
perty  data  for  Input  to  such  codes  was  made  with  the  Intention  of  minimizing 
the  complexity  of  the  analysis  while  remaining  on  firm  theoretical  ground. 

Basic  material  property  tests  which  would  supply  the  necessary  data  for  code 
Input  were  performed  using  a  simple  6061-T6  aluminum  tensile  specimen.  These 
tests  allowed  the  generation  of  "zero"  time  or  time  Independent  stress/straln/ 
temperature  curves  as  well  as  time  dependent  creep  strain  curves. 

Simulation  experiments  In  which  load  and  temperature  were  varied  as  a  function 
of  time  over  a  range  which  Included  stresses  beyond  yield  and  temperatures 
near  melt  (800°F)  were  performed  on  the  above  mentioned  tensile  specimens. 
Analyses  were  performed  using  MARC,  ANSYS  and  a  one-dimensional  code  "CREEPARHS" 
with  the  newly  generated  material  property  data.  The  analytical  results  were 
compared  with  the  results  obtained  from  the  simulation  experiments. 

2.2  CONCLUSIONS 

2.2.1  Mechanical  Material  Property  Data/Hodellnq 

Significant  time  dependent  creep  strain  (on  the  order  of  the  mechanical  strain 
or  greater)  can  occur  In  structures  where  stresses  exceed  yield  and  temperatures 
approach  melt,  even  In  the  1-10  second  engagement  time  regime.  It  thus  becomes 
necessary  to  account  for  the  time  dependent  and  time  Independent  effects.  This 


can  be  effectively  accomplished  by  developing  "zero"  time  or  time  independent 
stress/strain/temperature  curves  or  relations  augmented  by  time  dependent 
creep  strain  relations.  The  use  of  these  data  resulted  In  a  much  improved 
analytical/experimental  correlation  of  the  strain/time  response  in  simple 
uniaxial  simulation  experiments  as  compared  to  the  correlation  based  upon  the 
use  of  10  second  isochronous  or  standard  handbook  data.  This  is  because  the 
Isochronous  and  standard  data  have  built  In  ti. ’’e-in teg rated  creep  strain.  The 
new  data  is  particularly  useful  in  predicting  response  at  very  early  times 
(1-3  seconds)  and  will  In  fact  predict  much  higher  short  time  allowable 
stresses  for  a  one  dimensional  specimen.  An  investigation  into  the  effects  of 
varying  the  heat  up  time  or  the  time  required  for  the  specimen  to  reach 
test  temperature  prior  to  loading  showed  no  significant  trend  in  the  15-60 
seconds  heat  up  time  regime.  This  suggests  that  material  property  changes 
due  to  the  above  mentioned  heat  up  time  has  little  effect  upon  the  measured 
properties. 

2.2.2  Nonlinear  Structural  Response  Codes 

Theoretical  considerations  and  excellent  one  dimensional  experiment/analysis 
correlation  suggest  that  typical  general  purpose  codes  (MARC  and  ANSYS)  and 
specialized  codes  (CREEPARHS)  are  capable  of  analyzing  simple  and  complex 
structures  subject  to  rapid  heating  given  the  proper  material  data  input.  The 
codes  allow  separation  of  time  dependent  and  time  independent  material 
property  data  and  thus  are  well  suited  to  perform  combined  thermal  elastic 
creep  and  plasticity  analysis  as  required  in  the  rapid  heating  problems.  The 
general  purpose  codes  are  very  expensive  to  use  in  their  nonlinear  mode. 

Coding  errors  which  existed  in  the  nonlinear  temperature  dependent  subroutines 
in  MARC  and  ANSYS  have  been  corrected  but  it  is  suspected  that  these  codes  are 
still  not  error  free. 

2.3  RECOMMENDATIONS 

The  following  paragraphs  consist  of  a  list  of  specific  recommendations  (in 
order  of  suggested  performance)  for  future  work  followed  by  a  brief  discussion 
of  issues  inherent  in  the  recommendations.  These  recommendations  are  based 
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upon  results  from  present  and  previous  studies  concerned  with  the  ability  to 
predict  the  response  of  structures  subject  to  rapid  heating. 


Mechanical  Material  Properties 

•  Use  "new"  data  with  more  sophisticated  creep  laws  to  analyze  previously 
investigated  problems  such  as  a  centrally  heated,  bi-axial  loaded  plate. 

Two  Interrelated  matters  which  are  of  concern  in  determining  the  requirements 
for  material  property  data  are  constitutive  modeling  and  data  acquisition 
techniques.  The  modeling  of  the  data  should  be  sufficiently  sophisticated  to 
allow  for  the  many  combinations  of  possible  load/temperature/time  profiles  which 
the  structures  will  see  while  being  simple  enough  to  allow  for  straight¬ 
forward  raw  data  collection  as  well  as  ease  in  analysis.  The  thermally 
degraded  material  properties  can  have  a  significant  effect  on  the  predicted 
failure  times  and  failure  thresholds  of  even  simple  structures.  The  time 
dependent  or  creep  effects  must  be  Isolated  and  be  properly  accounted  for  in 
the  constitutive  modeling  and  analysis.  It  is  reconmended  that  time 
independent  "zero"  time,  piece-wise  linear  stress/strain  versus  temperature 
curves,  along  with  a  mechanical  equation  of  state  for  primary  and  tertiary 
creep  strain,  be  considered  as  the  basic  material  models.  In  addition,  various 
hardening  rules  (strain  hardening,  time  hardening,  etc.)  should  be  evaluated 
in  terms  of  their  predictive  capability. 

The  temperature  and  loading  environments  along  with  the  constitutive  models 
will  dictate  the  type  of  data  that  is  required.  It  then  becomes  necessary  to 
perform  experiments  which  will  provide  the  required  data.  In  the  present 
context,  it  is  necessary  to  bring  the  specimens  to  load  and/or  to  temperature 
in  a  very  short  period  of  time.  It  is  recommended  that  the  stress/strain/ 
temperature  data  be  obtained  by  bringing  the  specimens  to, and  stabilizing 
the  specimens  at  temperature  as  rapidly  as  possible  while  loading  the 
specimens  to  failure  in  a  time  period  which  is  small  compared  to  the  actual 
anticipated  time  to  failure  of  the  specimen.  Both  naw  anc  standard  data 
should  be  used  to  revisit  previously  analyzed  problems.  The  sensitivity  of 
tne  structure  to  a  wide  range  of  material  data  can  then  be  established. 
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<  Nonlinear  Code  Analysis 

•  Develop  advanced  closed  form  analysis  code  (CREEPARHS  II)  to  analyze 
more  complicated  problems 

-  Indeterminant  bar 

-  Circular  Plate 

•  Exercf se/val Idate  general  purpose  nonlinear  codes  on  above  mentioned 
probl ems. 

•  It  Is  essential  that  each  option  in  the  codes  be  exercised  and  vali¬ 
dated  prior  to  a  large  scale,  complicated  analysis.  It  Is  reconmended  that 
sample  problems  which  are  amenable  to  semi-closed  form  solutions  such  as  an 
indeterminant  bar  and  a  centrally  heated  circular  plate  be  used  to  verify  the 
operation  of  the  code  for  more  advanced  problems  using  the  "new"  data.  A  more 
advanced  version  of  the  CREEPARHS  code  could  serve  the  purpose  of  checking  the 
multidimensional  stress  calculational  ability  of  the  general  purpose  codes. 

Experiments 

•  Perform  appropriate  experiments  for  code/experiment  correlation. 

•  Simple  one  dimensional  simulation  tests  (load  and  temperature  vary 
with  time)  can  serve  as  the  vehicle  for  combined  code/materlal  property 
validation  for  more  refined  creep  and  hardening  laws. 

A  series  of  highly  instrumented  "simple"  experiments  which  demonstrate  the 
various  aspects  of  the  overall  rapid  heating  problem  should  be  performed. 

These  should  lend  themselves  to  straightforward  analysis,  and  can  be  limited 
to  elementary  structural  configurations. 

Failure  Laws 

•  Propose  and  evaluate  failure  laws  for  stress  and  buckling  type 
failures . 

Generate  Appropriate  Data  for  Alternate  Materials 

1.  References  for  and  discussion  of  the  need  for  "appropriate"  material 
properties . 

2.  Suggestions  for  types  of  simulation  tests. 

3.  Discussion  TTf  the  sensitivity  of  predicted  results  to  material 
properties/computer  code  modeling. 
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Short  Time  Heating  Analysis  Guide 

•  The  need  exists  for  a  manual  which  explains  the  various  aspects  of 
the  analysis  of  a  rapidly  heated  structure.  The  following  material  would 
be  Included  in  such  a  guide: 

1.  Summary  of  simplified  analysis  techniques  for  typical  structural 
configurations  Including;  bars,  plates,  shells. 

2.  Summary  of  relevant  experimental  data. 

3.  Appropriate  material  property  data  for  nonlinear  finite  element  code 
analysis. 


Section  3 

SURVEY  OF  ANALYTICAL  METIIODS/MODELS 

The  initial  phase  of  this  study  was  dedicated  to  developing  an  understanding 
of  the  methods  and  models  available  for  analyzing  the  nonlinear,  time 
dependent  behavior  of  metal:.  The  primary  consideration  was  to  evaluate  the 
suitability  of  existing  approaches, or  to  develop  an  alternative  approach  to 
material  modeling  which  could  be  easily  incorporated  into  existing  analysis 
codes.  The  material  model  must  have  a  firm  theoretical  base  but  be  simple 
so  as  to  minimize  the  amount  of  material  property  data  that  must  be  aathered. 

3.1  THEORETICAL  CONSIDERATIONS 

In  order  to  understand  the  behavior  of  metals,  for  example  6061 -Aluminum,  which 
are  rapidly  heated,  one  must  realize  that  time  dependent  plastic  flow  or  creep 
can  occur  even  in  very  short  periods  of  time  if  the  temperatures  and  stresses 
in  the  material  are  sufficiently  high.  Any  attempt  to  use  constitutive  models 
which  do  not  account  for  the  time  dependent  strain  explicitly  or  account  for 
some  time  integrated  effect  of  the  time  dependent  strain  (Isochronous  stress/ 
strain  curves)  is  likely  to  Introduce  major  errors  if  the  duration  of  heating 
or  loading  is  short  as  compared  to  the  time  over  which  the  data  are  taken.  In 
the  case  where  materials  are  heated  to  temperatures  approaching  melt  in 
less  than  10  seconds,  the  "appropriate"  data  must  account  for  the  strain  as 
a  function  of  temperature,  stress  and  time  for  times  less  than  10  seconds. 

Historically,  constitutive  models  for  homogeneous  structural  materials  were 
developed  In  a  manner  which  allowed  separation  of  the  material  behavior  into 
linear  (elastic)  and  nonlinear  (plastic)  time  independent  behavior  and  non¬ 
linear  time  dependent  behavior  (creep).  An  excellent  sunmary  discussion  of 
basic  concepts  of  plasticity  and  creep  theory  can  be  found  in  Reference  4. 

Table  1  Illustrates  the  various  regimes  of  material  behavior  as  is  used  in 
most  analysis  codes.  As  pointed  out  In  reference  4,  there  have  been  a  number 
of  attempts  to  Improve  on  the  classical  phenomenological  or  equation  of  state 
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method  of  describing  creep  strain.  These  theories  attempt,  for  example,  to 
account  for  certain  recovery  effects  found  upon  specimen  unloading.  Host 
recently,  the  so  called  “unified"  theories  based  upon  thermodynamic  consider¬ 
ations  have  surfaced.  These  theories  attempt  to  treat  the  inelastic  time 
dependent  and  time  independent  strain  as  a  single  phenomenon. 

The  point  to  keep  in  mind  is  that  the  particular  constitutive  model  chosen 
should  explain  the  observed  phenomenon  in  the  simplest  possible  manner.  In 
a  rapidly  heated  structure,  the  stress  redistribution  and  subsequent  failure 
arise  primarily  from  the  thermally  degraded  elastic  properties,  the  thermally 
activated  plastic  flow  and  the  early  time  creep.  These  phenomenon  are  pre¬ 
cisely  those  which  can  be  modeled  by  current  general  purpose  nonlinear  struc¬ 
tural  analysis  codes.  In  view  of  this,  the  material  behavior  will  be  modeled 
as  elastic-plastic-creep.  Selected  flow  ana  hardening  rules  will  be  used  to 
describe  the  plasticity  and  a  mechanical  equation  of  state  will  be  used  to 
model  the  creep. 

In  principle,  one  can  talk  about  time  independent  strain  (elastic-plastic)  and 
time  dependent  strain  (creep)  but  in  practice  ons  must  be  able  to  perform 
material  property  tests  in  a  manner  which  will  allow  separation  of  these 
quantities.  Naturally,  it  takes  a  finite  time  to  heat  a  specimen  to  temper¬ 
ature  and  aoply  the  load.  The  standard  test  to  obtain  stress/strain  data  as 
a  function  of  temperature  is  to  heat  the  specimen  to  temperature  in  a  minute 
or  longer  and  then  to  load  the  specimen  to  failure  at  a  low  strain  rate  of 
O.Wminut.s.  Obviously,  a  large  component  of  time  dependent  strain  can 
accumulate  during  the  test.  Theoretically,  one  ccula  obtain  "zero11  time  or 
time  independent  elastic/plastic  temperature  dependent  stress/strain  curves 
by  either  heating  a  specimen  to  temperature  and  applying  t.he  load  in  a 
very  short  period  of  time  or  by  applying  the  load  to  a  specimen  and  heating 
to  temperature  in  a  very  short  period  of  time.  The  methods  ars  essentially 
equivalent  in  the  case  where  the  heating  of  hr?  specimen  does  not  cause 
metallurgical  changes  in  the  time  it  takes  to  hea..  Tne  equivalence  of  the 
two  methods  has  been  a  subject  of  debate.  A  simple  thought  experiment  ca 
help  to  iliustiu^e  the  problems  in  laving  tv:o  sets  of  con.rlicting  "zero  t*m» 
data  as  generated  by  the  above  mentioned  approaches.  Consider  a  flight 
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structure  operating  at  some  nominal  temperature  and  subject  to  a  quasi-static 
set  of  aerodynamic  loads.  At  some  point  in  time,  the  structure  is  irradiated 
by  a  heat  source  and  portions  of  the  structure  are  rapidly  heated.  At  that 
instant,  one  might  argue  that  the  appropriate  data  might  be  that  generated  by 
applying  a  load  to  a  specimen  and  then  rapidly  heating  it.  Back  at  the 
structure,  in  the  next  instant,  the  local  stress  field  has  changed  due  to 
thermally  degraded  material  properties.  One  might  then  argue  that  the  appro¬ 
priate  data  might  be  that  generated  in  a  test  where  temperature  is  held 
constant  and  load  is  rapidly  changed  .  From  then  on,  the  stress  and  temperature 
are  in  a  constant  state  of  change.  It  is  obvious  that  neither  method  of 
obtaining  data  can  be  used  to  exactly  model  the  phenomenon  unless  the  differences 
between  the  data  are  negligible.  It  is  suggested  that  the  primary  difference 
in  the  data  will  arise  from  the  accumulated  time  dependent  strain  and  that 
either  method  will  produce  equivalent  "zero"  time  curves  as  long  as  the 
loading  or  heating  times  are  short  compared  to  the  engagement  time. 

3.1.1  Determination  of  "Zero"  Time  Temperature  Dependent  Stress-Strain  Curves 
A  procedure  has  been  developed  which  allows  the  determination  of  the  "zero" 
time  strain  from  stress-strain  tests  performed  at  various  load  rates.  Consider 
a  uniaxial  specimen  as  shown  in  Figure  2.  Figure  3  shows  the  load  and  temper¬ 
ature  history  for  the  specimen. 
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Figure  2.  Uniaxial  Test  Specimen 
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Figure  3.  Material  Property  Test  for  Constant  Temperature  Rapidly  Loaded 
Specimen 


The  specimen  is  uniformly  heated  in  20  seconds  or  less  to  test  temperature 

dP 

T-|  and  then  loaded  rapidly  at  a  rate  P-j  =  -j-  to  final  load  PQ.  Care  must  be 
taken  so  as  to  load  the  specimen  rapidly  enough  to  prevent  time  dependent 
strain  from  accumulating  but  not  so  rapidly  as  to  produce  dynamic  effects. 

The  strain  is  measured  for  each  load  P,  0  ^  P  _<  PQ  for  fixed  load  rate  P-|  and 
temperature  T^ .  A  series  of  tests  are  performed  in  which  strain  is  measured 
for  various  load  rates  and  for  particular  loads  (stresses)  and  temperatures. 
The  strain  is  then  plotted  as  a  function  of  the  time  to  reach  a  particular 
load  as  shown  in  Figure  4.  These  points  are  fit  by  a  smooth  curve  and  extra¬ 
polated  back  to  time  zero.  This  then  gives  the  time  independent  elastic- 
plastic  strain  for  a  given  stress  and  temperature.  A  series  of  tests  at 
various  appropriate  load  rates  and  temperatures  would  provide  data  that  can 
be  used  to  generate  "zero"  time,  temperature-dependent  stress-strain  curves. 


At 

TIME  TO  LOAD  (SEC) 

Figure  4.  Determination  of  "Zero"  Time  Elastic-Plastic  Strain 


15 


3.2  NONLINEAR  STRUCTURAL  ANALYSIS  CODES 

Many  general  purpose  and  special  purpose  codes  exist  which  have  the  capability 
of  solving  elastic-plastic  creep  oroblems.  The  solution  algorithms,  however, 
differ  In  many  respects  both  from  a  theoretical  and  operational  point  of  view. 
An  excellent  description  of  the  various  aspects  of  the  solution  techniques 
can  be  found  in  Reference  5.  The  thermal  elastic-plastic  algorithms  are 
fairly  standard  and  are  flexible  enough  to  allow  for  various  types  of  material 
behavior  beyond  yield  and  for  reversed  loading.  The  solution  methods  for 
handling  the  creep  strain  are  varied  and  can  greatly  influence  the  accuracy 
of  the  results  and  the  computation  time  of  the  analysis.  The  original  solution 
methods  (and  those  used  most  widely  today)  were  based  upon  the  method  used 
for  elastic  thermal  stress  analysis.  In  that  approach,  a  set  of  "initial" 
strains  depending  on  temperature  and  thermal  coefficient  of  expansion  are 
calculated  and  con\erted  to  "initial"  loads  through  the  use  of  elastic  and 
geometric  properties.  The  problem  is  linear  since  none  of  the  "initial" 
quantities  depended  on  the  deformation  or  stress  level  in  the  structure. 
Similarly,  "initial"  creep  strains  can  be  calculated.  These  strains  in  general 
depend  on  the  stress  level  of  the  structure  and  thus  can  only  serve  as 
approximate  values.  Many  alterations  to  the  basic  "initial"  strain  procedure 
have  been  made  to  account  for  the  approximate  nature  of  the  solution,  including 
algorithms  which  select  time  increments  so  as  to  insure  that  stresses  calcu¬ 
lated  at  the  end  of  an  increment  are  compatible  with  creep  strain  calculated 
at  the  beginning  of  the  increment. 

3.2.1  General  Purpose  Finite  Element  Codes  (MARC,  ANSYS) 

Two  general  purpose  codes  were  selected  to  demonstrate  and  evaluate  the 
ability  of  standard  nonlinear  codes  to  solve  the  rapid  heating  problems.  These 
codes  have  extensive  element  libraries  and  are  capable  of  solving  nonlinear 
problems  including  plasticity,  creep,  large  deflections  and  buckling.  The 
codes  use  essentially  the  same  approach  in  solving  the  plasticity  and  creep 
problems  although  MARL  allows  for  user  written  subroutines  to  calculate  creep 
and  temperature  dependent  stress-strain  behavior.  A  more  detailed  description 
of  the  codes  can  be  found  in  References  3  and  5. 
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3.2.2  Special  Purpose  Analysis  Code  (CREEPARHS) 


A  special  purpose  one-dimensional  direct  Integration  code  "CREEPARHS"  (creep- 
elastic-plastic-analysis  of  rapidly  heated  specimens)  was  written  for  the 
expressed  and  Important  purpose  of  evaluating  the  proposed  material  models 
In  an  accurate  and  cost  efficient  manner.  The  code  is  capable  of  analyzing 
the  response  of  an  idealized  one-dimensional  homogenous  cylindrical  bar 
subjected  to  load/ temperature  histories  as  shown  in  Figure  5. 


Figure  5.  Load/Temperature  Input  Histories  for  Simulation  Tests 

These  temperature  and  load  histories  were  selected  as  being  representative  of 
those  seen  in  more  complex  rapidly  heated  structures  and  to  demonstrate  the 
wide  range  of  material  behavior  that  can  result  from  varying  the  input  load 
and  temperature  levels  and  durations. 

The  code  uses  time  independent  multilinear  stress-strain-temperature  curves 
and  interpolates  for  temperatures  between  those  Inputs.  In  addition,  it  cal¬ 
culates  creep  strain  as  a  function  of  time  using  a  specific  creep  strain  versus 
stress,  time  and  temperature  relations,  and  a  time  hardening  creep  rule. 

Finally,  thermal  strain  as  a  function  of  temperature  Is  calculated.  The 
output  consists  of  temperature,  stress,  and  mechanical,  thermal,  creep  and 
total  strain  at  each  point  In  time  as  well  as  a  listing  of  the  Input  data  and 
initial  conditions.  A  more  detailed  description  of  the  code  including  the 
governing  equations  is  given  in  Section  7.  A  code  listing  can  be  found  In 
Appendix  2.  It  should  be  noted  that  the  input  to  this  code  Is  identical  to 
that  required  by  most  general  ourpose  codes.  This  allows  a  direct  check  on  the 
various  procedures  used  by  the  codes  to  solve  the  nonlinear  problems. 
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Section  4 

EXPERIMENT/TEST  SPECIMEN  DESIGN  AND  DATA  REDUCTION 

Simple  uniaxial  tensile  specimens  were  used  for  all  of  the  material  property 
and  simulation  tests.  The  material  property  tests  consisted  of  stress/straln 
tests  as  a  function  of  temperature  and  short  time  creep  tests  as  a  function 
of  stress  and  temperature.  The  experiment  and  specimen  designs  were  based 
upon  ASTM  El 50-64,  standard  recommended  practices  for  "conducting  croep  and 
tension  tests  of  metallic  materials  under  conditions  of  rapid  heating  and 
short  time". 

Subsection  4.1  through  4.3  discuss  the  tests  and  test  specimens.  Subsection 
4.4  describes  the  process  used  to  convert  the  raw  data  to  a  form  suitable  for 
input  to  the  analysis  codes.  This  was  a  major  effort  and  should  be  automated 
In  any  future  work. 

4.1  SPECIMEN  DESIGN 

The  tensile  specimens  were  machined  from  a  sheet  of  6061-T6  Aluminum,  The 
final  specimen  dimensions  (Figure  6)  were  determined  experimental ly  by  varying 
the  specimen  length  and  cross  section  In  order  to  obtain  the  most  uniform 
temperature  field  In  the  gauge  length. 


4.50"  (11.43  cm) 
2~ 7 5"  (6.90  cm) 

_ _ _ L__ _ 


1 

I  0,200" 


1 


( .  71  cm)D1*.  0.37 


( , 90  cm)  D1«  . 


Figure  6.  Tensile  Specimen  Design 
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4.2  TEST  EQUIPMENT 

The  following  Itema  wort  used  In  tht  rapid  ht&tlng/loadlng  teats: 

1.  MTS  model  810  Servo- hydraulic  teat  machine  with  a  100,000  pound 
load  frame  end  a  5,000  pound  load  cell. 

2.  MTS  model  632.558*01  elevated  temperature  extenaometer  with  quartz 
arma  and  a  one-inch  gauge  length. 

3.  Reaearch  Inc.  THERMAC  "A"  Power  controller  with  DATA-TRAC  Program 
Control.  100  AMP,  440  volt  capacity. 

4.  Jefferlea  Tranaformer,  460  volt  primary,  16  volt  secondary,  15  KVA 
capacity. 

4.3  TEST  TECHNIQUE 

Direct  resistance  heating  was  used  In  order  to  obtain  the  rapid  heating  rates 
required  In  the  testa.  A  16  KVA  transformer  coupled  with  a  Reaearch  Inc. 

Power  Controller  equipped  with  a  date-trac  programmer  was  used  to  control 
heating  rates  and  maximum  temperatures,  The  specimens  were  Instrumented  with 
4  ci.rome-Alumal  thermocouples  attached  by  the  use  of  copper  clamps  made  from 
thin  copper  shoeti.  One  thermocouple  was  placed  on  each  side  of  the  one-inch 
gauge  length  to  monitor  th*  temperature  variation.  The  remaining  two  thermo¬ 
couples  were  used  to  monitor  and  control  the  temperature  within  the  gauge  length. 
The  maximum  heating  rate,  maximum  temperature,  loading  rate  and  maximum  load 
were  all  preset  on  the  test  equipment  prior  to  the  start  of  each  test.  A 
diagram  of  the  test  equipment  Is  shown  In  Figure  7.  Figures  8  and  9  are  photo¬ 
graphs  showing  the  entire  test  system  and  a  closeup  of  the  test  speclmen/geuge 
configuration,  Many  of  the  tests  required  heating  prior  to  loading.  For 
these  tests,  a  zero  offset  technique  was  used  to  offset  the  strain  trace  by 
a  known  amount,  This  allowed  the  use  of  a  higher  magnification  factor  for 
the  strain  versus  time  oscillograph  trace  during  the  loading  phase  of  the  test, 

The  following  Is  a  list  of  the  range  of  loading  and  heating  parameters  which 
were  used  In  the  tests: 

e  load  rate  (?.6  -  HO  K5I/5ec) 
e  stress  (0-45  K51) 
e  Mealing  rate  (1b  -  4Q*f/5»ec) 

e  Temp  variation  across  gauge  length  (O0"f  at  H00*F) 
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4.4  DATA  OUTPUT  AND  REDUCTION 

The  strain  versus  time  as  well  as  the  temperature  and  load  histories  were 
mechanically  recorded  on  a  standard  continuously  running  oscillograph  recorder. 
In  addition,  the  stress-strain  curves  were  plotted  in  real  test  time  using  an 
X-Y  plotter.  An  example  of  the  o-graph  output  is  shown  in  Figure  10. 

A  major  effort  in  the  program  was  to  convert  the  recorded  data  into  a  form 
which  was  amenable  to  analysis.  The  data  reduction  task  was  complicated  by 
the  fact  that  the  strain  scales  and  chart  recording  speeds  were  varied  up  to 
three  times  during  a  test  and  from  test  to  test.  The  scale  changing  was 
necessary  because  the  magnitude  of  the  strain  data  varied  greatly  both  during 
the  tests  and  from  test  to  test.  The  strain  versus  time  was  digitized  using 
a  Tektronix  digitizer  and  then  stored  on  magnetic  tape. 

A  computer  program  was  written  which  smoothed  and  provided  curve  fitting  for 
the  raw  data  and  then  calculated  the  parameters  necessary  for  analysis  such 
as  the  elastic  moduli,  yield  stress,  and  stress  versus  plastic  strain.  In 
particular,  the  stress-strain  data  was  fit  to  a  10th  order  polynomial  and 
then  reduced  to  a  piecewise  linear  curve  with  N  *ej  lents  (Figure  11).  The 
N  +  1  breakpoints  were  selected  by  calculating  the  N  points  on  the  curve  with 
the  largest  changes  in  slope.  The  idea  was  to  best  approximate  the  work 
hardening  effects  In  the  stress  strain  behavior. 

The  creep  data  was  fit  by  various  formulas  depending  upon  the  intended  use. 


The  general  behavior  of  the  creep  data  Is  illustrated  In  Figure  12.  One  sees 
the  classical  primary,  secondary  and  in  a  few  cases  tertiary  steps.  In  most 
cases,  the  time  scale  of  the  secondary  stage  was  very  long  compared  to  1-10 
second  heating  time  which  Is  of  current  interest  and  thus  the  data  was  fit 
to  a  fourth  order  polynomial  using  a  least  squares  technique  In  the  primary 
stage  and  was  assumed  linear  in  the  secondary  stage.  The  curve  fit 
representations  of  the  creep  data  are  shown  in  the  next  section  (Figures 
20  -  24). 
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Figure  12.  General  Behavior  of  Creep  Strain  Versus  Time 


Section  5 

MATERIAL  PROPERTY  TESTS 

This  section  Includes  a  description  of  the  tests  and  a  surrmary  of  the  test 
results  followed  by  a  complete  test  matrix.  The  emphasis  In  this  program  was 
not  to  obtain  statistically  significant  data,  but  Instead  to  perform  tests 
which  would  Illustrate  basic  material  behavior.  Reruns  of  individual  tests 
were  performed  only  In  those  cases  where  the  data  did  not  follow  anticipated 
trends.  In  all  cases  where  the  data  are  used  In  an  analysis,  the  error  in 
the  data  Is  considered  small  compared  to  the  magnitude  of  the  effect  being 
measured. 

The  following  types  of  tests  were  performed  to  obtain  basic  material  property 
data.  A  sketch  of  the  load  and  temperature  history  for  the  test  is  included 
In  each  description. 

5.1  TEMPERATURE  PREHEAT  SENSITIVITY  TEST-' 

These  tests  were  perform^ _  t?  gain  an  insight  into  the  effect  of  varying  the 
preheat  time  on  the  creep  maracterii ,1cs  of  the  material.  In  particular, 
the  time  to  reach  800°F  i427°C)  was  varied  from  2G-C0  seconds.  The  lead  was 
then  applied  at  a  rate  of  40  KSI/second  (275  MPa/second)  to  a  final  stress  of 
4  KSI  (27.5  MPa).  The  strain  as  a  function  of  time  was  measured  at  constant 
stress.  The  load  and  temperature  histories  are  shown  In  Figure  13. 

Tne  results  of  these  tests  are  shown  In  Figure  14.  The  cross  hatched  region 
In  the  lower  curves  indicates  the  area  where  the  strain  versus  time  data  lies. 
There  was  no  trend  noticed  In  the  behavior  of  the  strain  versus  time  as  the 
haat-uo  time  was  varied.  This  is  not  to  say  that  no  effect  exists.  The 
upper  curve  was  a  similar  single  test  performed  at  40  KSI/second  (275  MPa/ 
second)  with  a  20  second  heat  up  time.  It  Is  suggested  that  a  variation  of 
a  principal  quantity  such  as  stress  Is  far  more  significant  In  determining 
the  strain-time  behavior  than  a  variation  In  heat  up  time.  A  test  matrix  is 
given  ■’n  Table  2. 
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Figure  13.  Stress/Tempereture  History  -  Preheat  Sensitivity  Tests. 


Table  2 

Preheat  Sensitivity  Test  Matrix 


Test 

Type 

Test  No. 

Heatup 

Rate 
( °F/Sec ) 

Final 

Stress 

(KSI/MPa) 

CREEP 

CT1 -80 

40 

4/27,5 

CREEP 

CT2-80 

32 

4/27.5 

CREEP 

CT3-80 

16 

4/27.5 

CREEP 

CT4-80 

8.9 

4/27.5 

5.2  STRAIN  VERSUS  LOAD  RATE  TESTS 

These  tests  were  useu  to  demonstrate  the  ability  to  extrapolate  the  measured 
strain  to  zero  time  so  as  to  generate  the  "zero"  time  or  time  independent 
s tress- strain- temperature  curves. 

The  specimens  were  tested  at  000°F  (260°C)  end  800°F  (427°C)  at  stresses  well 
above  yield  so  as  to  be  able  to  bound  the  time  dependent  effects  due  to 
temperature  anc  stress.  The  general  procedure  was  to  heat  the  specimens  to 
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Figure  14.  Strain  vs  Heatup  Time  -  Temp  =  800°F  (427°C) 


temperature  Tq  In  20  seconds  and  then  apply  the  load  at  a  constant  rate 
until  the  final  load  level  was  reached.  The  load  rate  was  then  changed 
to  and  the  test  repeated.  These  histories  are  shown  In  Figure  15.  The 
specific  tests  performed  are  listed  in  Table  3. 


Figure  15.  Stress/Temperature  History  -  Load  Rate  Sensitivity  Tests. 


Table  3 

Load  Rate  Sensitivity  Test  Matrix 


Test 

Type 

Test  No. 

Test 

Temperature 

°F/°C 

Stress  Rate 
(KSI/MP,j/Sec ) 

Final 

Stress 

(KSI/MPa) 

Cr< 

»ep 

CS1-50 

500/260 

80/550 

28.7/198 

CS2-50 

500/260 

20/137 

28.7/198 

CS3-50 

500/260 

10/69 

28.7/198 

CS1-80 

800/427 

40/275 

4/27.5 

CS2-80 

800/427 

20/137 

4/27.5 

CS3-80 

800/427 

10/69 

4/27.5 

CS4-80 

800/427 

2.5/14 

4/27.5 

CS5-80 

800/427 

40/275 

6/41 .3 

CS6-80 

800/427 

20/137 

6/41  .3 

CS7-80 

800/427 

10/69 

6/41.3 

<  The  strain  at  500°F  (260°C)  is  plotted  as  a  function  of  the  time  to  reach  the 
final  stress  (28.7  KSI/198  MPa)  In  Figure  16.  A  smooth  line  is  drawn  through 
the  data  points  and  extrapolated  to  time  t  =  0.  The  Intercept  strain  repre¬ 
sents  the  time  Independent  Elastic  +  Plastic  strain  at  the  given  stress  and 
temperature.  It  is  seen  that  a  reasonable  approximation  to  the  time  independent 
strain  (<15%  error)  can  be  obtained  by  loading  the  specimen  to  final  stress  in 
one  second  or  less.  On  the  other  hand,  a  200%  error  would  result  if  the 
specimen  were  loaded  in  three  seconds.  Similar  results  are  shown  in 
Figure  17  for  specimens  tested  at  800°F  (427°C)  at  two  stress  levels.  It  is 
seen  tnat  for  a  given  acceptable  level  of  error,  the  time  tc  final  load 
increases  dramatically  as  the  stress  increases  beyond  yield.  In  this  case,  a 
15%  error  or  less  is  guaranteed  by  loading  to  6  KSI/41  MPa  in  0.15  seconds 
or  less.  Theoretically,  one  could  perform  a  whole  series  of  tests  where  the 
final  load,  load  rate  and  temperature  were  varied.  Extrapolation  of  the 
results  would  lead  to  a  set  of  "zero"  time  stress/strain-temperature  curves. 

In  practice  one  could  decide  upon  an  acceptable  upper  error  bound  and  select 
a  single  load  rate  which  would  be  used  In  all  of  the  tests.  In  this  study,  a 
load  rate  of  40  KSI/sec  (276  MPa/sec)  was  chosen  based  upon  a  maximum  15% 
error  in  the  strain  at  the  worst  case  teinperature/stress  combination 
(800°F/6  KSI)(427°F/41  MPa).  This  will  assure  a  minimum  error  in  the 
majority  of  the  stress/strain/temperature  data. 

5.3  STRESS-STRAIN-TEMPERATURE  TESTS 

The  tests  that  were  performed  in  this  phase  were  essentially  standard  elevated 
temperature  stress-strain  tests  with  the  exception  that  the  load  was  applied 
at  a  rate  (40  KSI/sec)  (276  MPa/sec)  which  was  rapid  enough  to  prevent  the 
accumulation  of  time  dependent  strain  (see  Section  5.2).  Fig.re  18  illustrates 
typical  load  temperature  histories.  Table  4  contains  the  complete  test 
matrix. 

The  resulting  "zero"  time  stress-strain  curves  are  plotted  as  a  function  of 
temperature  in  Figure  19.  It  is  important  to  note  that  a  relatively  small 
increment  of  stress  beyond  yield  at  temperatures  above  500°F  (260°C)  leads 
to  an  increasingly  large  increment  in  plastic  strain.  An  important  comparison 
between  these  curves  and  the  so  called  handbook  and  isochronous  curves  is 
presented  in  the  following  section. 
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Streia*5tri1n  Tamparatura  Tatt  Matrix 


T«it  Typa 

Tilt  No. 

Teat 

Tamparatura 

(Vc) 

Final  Strati 
(KSl/MPa) 

Tima  To 
Final  Strati 
(Stcoridt) 

Stroii/Straln 

SS1-RT 

71/21.7 

42/289 

1 .06 

f 54 -40 

400/204 

33/227 

0.025 

556-60 

500/260 

20,7/198 

0,7176 

556-60 

600/315 

20/138 

0,6 

557-70 

700/371 

11/76 

0.275 

SSO-76 

760/399 

0/5,6 

0,20 

559-00 

000/42/ 

6/4,3 

0,16 

t».4  IRUP  TLSTS 

A  tar1at  of  ttandard  craap  taitt  wara  parformad  In  which  tht  load  and  tamper- 
•  ttjra  tail  lavalt  war*  ulactad  to  at  to  allow  Invattlyatlon  of  tht  cftop 
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response  of  the  material  over  a  wide  range  of  loading  conditions.  In 
particular,  the  temperature  levels  varied  from  500-800°F  (260-427°C)  and 
the  load  levels  were  selected  so  that  at  least  one  test  was  performed  in  the 
following  stress-strain  regions: 

•  elastic 

•  yield 

•  post-yield 

The  specimens  were  heated  to  temperature  in  20  seconds  or  less  and  the  load 
was  applied  at  the  rate  of  40  KSI/second  (276  MPa/second).  The  temperature/ 
load  history  Is  the  same  as  that  seen  in  Figure  13.  The  complete  creep  test 
matrix  is  shown  in  Table  5. 


Table  5 

Creep  Strain  Test  Matrix 


Test  Type 

Test  No. 

Test  Temperature 

Test  Stress  (KSI/MPa) 

Creep 

Cl -50 

500°F/260°C 

13/90  ! 

C2-50 

20/138 

C3-50 

25/172 

C4-50 

27/186 

Cl -60 

600°F/316*C 

8.1/56 

C2-60 

11.4/79 

C3-60 

17.8/123 

C4-60 

19.5/134 

Cl -70 

700°F/371 °C 

5.5/38 

C2-70 

6.8/47 

C3-70 

10.3/71 

C4-70 

■ 

11.3/78 

Cl -75 

750°F/399°C 

4.5/31 

C2-75 

5.8/40 

C3-75 

7.8/54 

C4-75 

8.4/58 

Cl -80 

800°F/427°C 

1.6/11 

C2-80 

2.4/17 

C3-80 

3.2/22 

' 

C4-80 

4.0/28 
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The  creep  strain  versus  time  curves  were  obtained  by  subtracting  the  strain  at 
load  and  then  curve  fitting  the  data  as  described  in  Section  4.4.  The 
resulting  curves  are  shown  In  Figures  20  -  24.  Note  that  creep  strain  increased 
dramatically  once  a  certain  stress  level  is  reached.  This  stress  level  is 
generally  In  the  center  of  the  knee  of  the  stress-strain  curve  as  can  be  seen 
from  Figure  19.  One  then  must  consider  the  contribution  from  creep  at  stress 
levels  above  yield. 

The  raw  data  from  the  creep  tests  (which  includes  elastic  and  plastic  strain) 
can  also  be  used  to  generate  so-called  Isochronous  or  constant  time  stress- 
strain  curves.  It  Is  precisely  these  isochronous  curves  that  have  been  used 
in  attempts  (Reference  3)  to  solve  the  rapid  heating  problems.  Unfortunately, 
the  best  available  Isochronous  data  are  10  second  curves.  It  is  evident 
from  the  creep  data  that  significant  creep  can  occur  during  the  1-10  second 
time  period,  thus  the  10  second  data  are  not  appropriate  for  the  analysis 

The  isochronous  curves  are  generated  by  selecting  a  certain  time  and  then 
plotting  the  stress  versus  strain  as  a  function  of  temperature.  An  accurate 
isochronous  curve  thus  requires  data  from  the  large  number  of  creep  tests. 
Examples  of  Isochronous  curves  which  have  been  generated  from  the  creep  data 
In  Figures  20-24  are  shown  in  Figures  25-29.  The  Important  thing  to  notice 
in  these  curves  Is  that  for  a  particular  value  of  strain,  the  corresponding 
stress  decreases  dramatically  as  one  moves  from  the  1  second  to  the  10  second 
curve.  If  the  10  second  data  were  used  in  the  analysis  of  a  simple  uniaxial 
specimen,  one  would  predict  failure  at  a  stress  level  which  could  be  signifi¬ 
cantly  lower  than  the  actual  failure  stress  level.  This  suggests,  for  example, 
that  the  predicted  failure  times  In  more  complicated  structures  might  be  much 
shorter  than  the  actual  failure  times. 

5.5  FREE  THERMAL  EXPANSION  TESTS 

In  each  of  the  previous  tests,  the  thermal  expansion  strain  was  measured  after 
the  specimen  was  stabilized  at  temperature.  The  coefficient  of  thermal 
expansion  was  calculated  by  dividing  the  thermal  strain  by  the  temperature. 

These  data  were  used  as  Input  to  the  analysis  codes. 
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Figure  21.  CREEP  STRAIN  VS  TINE  -  TERP  =  600°F  (316°C) 
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Figure  26.  ISOCHRONOUS  STRESS-STRAIN  CURVE  -  TEMP  =  600°F  (316°C) 
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•4  6.6  STRESS-STRAIN  CURVE  COMPARISON 

The  major  fault  In  using  either  10  second  Isochronous  or  standard*  (Reference 
6)  stress-strain  curves  It  that  a  substantial  amount  of  tlma  Integrated  creep 
strain  It  automatically  built  In.  This  can  ba  aaslly  saan  by  comparing  the 
stress-strain  temperature  curvet  at  shown  In  Figure  30.  The  first  thing  to 
notice  It  that  the  10  second  Isochronous  curves  do  not  differ  significantly  from 
the  standard  curvet,  especially  at  temperatures  above  500*F  (260*C).  However, 
these  curves  are  substantially  different  than  the  "zero"  time  curves.  For 
example,  at  600*K  (316*C)  there  It  almost  a  factor  of  two  difference  In  the 
maximum  stress.  One  must  not  forget  to  add  the  appropriate  creep  strain  to  the 
"zero"  time  strain  to  get  a  precise  comparison  between  stress-strain  behavior, 
but  this  correction  will  become  lets  significant  at  the  time  at  twnperature  of  the 
specimen  It  reduced.  In  the  limits  1 .e. ,  at  the  specimen  It  Instantaneously 
brought  to  temperature  and  load, the  creep  effect  would  vanish  and  the  "zero" 
time  data  must  be  used. 


Standard  (handbook)  curves  ere  typically  Generated  under  low  load  rate 
(0,b%  stra1n/m1mite)/lony  soak  time  (30  minute*..)  conditions. 
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Figure  30.  STRESS/STRAIN  CURVES:  "ZERO  TIME"/ 10  SEC  ISOCHRONOUS/STANDARD 


Section  6 

SIMULATION  EXPERIMENTS 


A  rapidly  heated  complex  structure  would  experience  a  continuous  redistribution 
of  temperature  and  stress  due  to  heat  flow  and  material  property  thermal 
degradation.  In  order  to  have  confidence  in  the  ability  of  an  analytical 
procedure  to  predict  the  response  of  the  structure,  one  could  perform  experi¬ 
ments  which  simulate  the  important  aspects  of  the  structural  response  and  then 
compare  the  results  with  the  analytical  prediction.  To  this  end,  a  series 
of  uniaxial  simulation  experiments  were  performed  in  which  the  load  and  temper¬ 
ature  histories  were  selected  so  as  to  cover  the  spectrum  of  the  postulated  load/ 
temperature  histories  in  the  actual  rapidly  heated  structures.  The  test  times 
were  slightly  longer  (<25  seconds)  than  the  suggested  1-10  second  rapid  heating 
times.  These  times  were  chosen  so  as  to  test  the  suitability  of  the  creep 
relation  for  longer  times  and  so  as  to  be  able  to  reach  800°F  (427°C)  without 
overstressing  the  experimental  heating  system.  These  experiments  served  three 
important  functions.  The  first  was  to  demonstrate  the  ability  to  perform  the 
simulation  experiments  in  a  precisely  controlled  manner.  The  second  was  to 
provide  an  experimental  base  which  could  be  used  in  comparison  with  analytical 
results.  The  third  was  to  allow  the  investigation  of  material  behavior  as 
influenced  by  various  combinations  of  time  varying  stress  and  temperature. 

This  investigation  ultimately  would  allow  one  to  make  an  assessment  of  the 
adequacy  of  existing  and  proposed  material  models. 

The  most  general  form  of  the  load/ temperature  history  is  shown  in  Figure  31. 
Notice  that  these  forms  allow,  as  special  cases,  all  of  the  tests  performed 
in  Section  5.  Within  this  framework,  16  different  types  of  tests  were 
performed  in  which  the  following  conditions  were  enforced: 

•  Constant  load  and  temperature  rising  with  time 

•  Load  and  temperature  rising  with  time 

•  Increasing  load  followed  by  decreasing  load,  temperature  constant 

•  Increasing  load  followed  by  decreasing  load,  temperature  rising 
with  time. 


Time  Time 


Figure  31.  Load/Temperature  History  Simulation  Experiments 

The  following  subsections  will  describe  results  from  one  example  of  each  of 
the  above  categories.  The  balance  of  results  are  plotted  aginst  the  results 
from  the  CREEPARHS  analysis  and  are  included  in  Appendix  1.  A  complete  simulation 
test  matrix  Is  given  In  Table  6.  It  will  be  seen  that  any  combination  of 
thermal,  mechanical  and  creep  strain  components  can  be  the  primary  strain 
response  modes. 

6.1  SIMULATION  TEST  -  CONSTANT/LOAD  TEMPERATURE 

The  strain  versus  time  trace  along  with  the  stress-temperature-time  histories 
are  shown  in  Figure  32.  The  strain  is  composed  of  free  thermal  expansion, 
mechanical  and  creep  components.  From  the  stress-strain  curves  in  Figure  19, 
one  can  see  that  the  mechanical  strain  is  less  than  0.1%  until  the  temperature 
exceeds  750°F  (399°C).  We  thus  see  an  almost  linear  rise  in  strain  due  to 
thermal  expansion  until  18  seconds  when  the  stress  and  temperature  are  constant. 
This  temperature/stress  combination  causes  significant  creep  (Tertiary)  from  this 
point  on  and  In  fact  the  specimen  ultimately  failed.  In  order  to  predict 
the  total  behavior,  one  must  account  for  tertiary  creep  in  the  analysis. 

CONCLUSION:  Thermal  strain  and  tertiary  creep  are  dominant  modes. 

6.2  SIMULATION  TEST  -  LOAD/TEMPERATURE  RISING  WITH  TIME 

The  strain  versus  time  trace  Is  shown  in  Figure  33.  The  strain  is  again  composed 
of  thermal,  mechanical  and  creep  effects.  The  response  Is  essentially  linearly 
elastic  until  one  second  when  the  maximum  stress  of  20  KSI  (138  MPa)  is  acting. 
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Table  6.2.  Simulation  Parameters  -  Increasing  Load  /Increasing  Temperature 


Table  6.0.  Simulation  Parameters  -  Increas ing-Decreasing  Load/Constant  Temperature 


Table  6.4.  Simulation  Parameters  -  Increasing-Decreasing  Load  /Increasing  Temperature 


.STRESS 


From  Figure  19  one  can  see  that  the  response  is  elastic  for  20  KSI  (138  MPa) 
stress  until  the  temperature  exceeds  500°F  (260°C).  This  Is  seen  in  Figure 
33  as  the  thermal  and  elastic  strain  increase  almost  linearly  until  12  seconds 
have  elapsed.  As  the  temperature  exceeds  50Q°F  (260°C),  the  maximum  allowable 
stress  at  temperature  is  exceeded  and  the  specimen  rapidly  fails  due  to  plastic 
and  creep  strain. 

CONCLUSION:  Thermal  strain  and  plastic  strain  are  dominant  modes. 

6.3  SIMULATION  TEST  -  LOAD- INCREASING,  DECREASING/TEMPERATURE  CONSTANT 

This  is  the  first  example  of  a  case  where  unloading  occurs.  The  thermal  strain 
has  [due  to  constant  temperature  at  500°F  (260°C)]  been  subtracted  out  of  the 
strain-time  history  and  the  results  are  shown  in  Figure  34.  From  Figure  19,  one 
can  see  that  the  stress  is  always  in  the  elastic  range.  The  nonlinear 
behavior  seen  in  Figure  34  in  the  seven  to  ten  second  time  interval  must  be  due 
to  creep.  One  also  notices  that  unloading  causes  a  linear  reduction  in  strain. 
Finally,  a  permanent  component  of  strain  remains  after  the  load  has  dropped  to 
zero. 

CONCLUSION:  Elastic  and  creep  strain  are  dominant  modes. 

6.4  SIMULATION  TEST  -  LOAN- INCREASING,  CONSTANT,  DECREASING/TEMPERATURE 
INCREASING 

This  is  the  final  and  most  complicated  example.  Note  in  Figure  35  that  the 
maximum  stress  and  temperatures  are  the  same  as  in  the  previous  example  but 
the  strain  versus  time  history  is  very  different.  The  strain  is  composed  of 
elastic,  thermal  and  creep  components  in  this  case  The  creep  component 
plays  a  small  role  due  to  the  fact  the  temperature  does  not  rise  sufficiently 
until  the  end  of  the  test.  One  can  see  a  linear  increase  in  elastic  strain 
until  two  seconds  have  elapsed.  The  strain  increases  due  to  thermal  effects 
until  approximately  11  seconds  when  the  load  Is  reduced  to  zero.  The  strain 
reduces  linearly  and  only  the  thermal  and  a  small  creep  component  remain. 


CONCLUSION:  Elastic  and  thermal  strain  are  dominant  modes. 


TEMP 


Simulation  Test  3*  Load  (Increasing/Constant/Oecreas 
Tenp  Increasing 


Section  7 

ANALYSIS,  AND  EXPERIMENT  -  CORRELATION 

This  section  will  describe  the  results  f,h« t  were  obtelned  from  the  correlation 
analysla  and  compare  those  results  with  the  results  from  the  simulation  exper¬ 
iments.  In  particular,  four  examples  of  the  CREEPARHS  versus  test  data 
comparison  (one  for  each  general  category  of  test)  will  be  discussed.  Ihe 
difference  In  results  obtained  by  using  "zero"  time  data  and  10  second 
Isochronous  data  will  be  demonstrated,  A  comparison  of  the  results  obtained 
for  MARC,  AN5YS  and  CRECPAkMS  will  be  made  for  two  of  the  ebeva  mentioned 
problems. 

7.1  STRUCTURAL  CODE  "USABILITY"' 

It  It  appropriate  to  state  at  tha  outset  that  the  general  purposa  nonlinear 
codes  such  as  MARC  and  ANSYS  are  not  practical  for  use  in  Investigating  IN 
response  of  a  structure  to  meny  combinations  of  thermal  and  mechanical  applied 
load  states.  The  codes  are  expansive  to  use  (as  will  be  seen)  for  even  the 
simplest  of  problems,  In  many  cases,  the  dominant  physics  of  the  problem  are 
occurred  due  to  Urn  "black'  box"  approach  used  In  the  codes,  and  It  Is 
precisely  on  understanding  of  these  physics  that  offers  one  the  opportunity  to 
gain  an  understanding  of  the  effects  of  the  material  behavior  on  the  response 
of  a  structure. 

An  additional  and  ever  more  serious  problem  Is  that  general  purpose  coders  may 
not  be  error  fra*.  It  is  our  experience  that  the  H,?  version  of  MARC  has 
"bugs"  In  the  routlnas  which  allow  the  calculation  of  plastic  strain,  These 
"buys"  have  been  corrected  In  the  11,4  revision,  further,  we  found  that  the 
revision  3  update  37  version  of  AM'.Y'j  has  pro  hi  on-:  In  <  ulculatlng  the 
tempera  lure  dependent  silffness  for  the  first  •  IwimmI  In  in  y  model.  A  method 
of  circumventing  the  problem  allows  the  usn  of  this  version  In  the  'orient 
study,  for  them  reasons,  It  was  decided  to  develop  a  simple  one  dimensional 
code  GRLLPAKIlii  which  could  perforin  the  required  simulation  analysis,  Him 
majority  of  the  analysis  was  par  formed  using  this  rude. 
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7.2  MATERIAL  PROPERTY  MODEL  FOR  TENSILE  SPECIMEN 

Th«  tensile  specimen  used  In  tht  simulation  experiments  can  b«  characterized 
by  a  one-d  linens  lonel  stress  atate.  The  quantity  of  Interest  which  was  measured 
In  the  experiments  wai  the  extenslonil  strain  as  a  function  of  time.  The 
total  strain  can  be  decomposed  Into  elastic,  plastic,  thermal  and  creep 
components  at  shown  In  aquation  1. 

«  •  t£  (<>#t)  ♦  ip(T,t)  ♦  tTH(T)  ♦  *-cr(°.T,t)  (1) 

where 

?  elastic  strain 

rp  ?  plastic  (time  Independent  strain) 
tTH  tt  thermal  strain 
r_cr  s  creep  strain 

and  o#  T,  t  are  applied  stress,  temperature,  and  time. 

The  elastic  component  can  be  rewritten  as 

*f  ’  rfl)  12) 

where  t'.(T)  U  the  elastic  temperature  dependent  Young's  modulus.  The 
elastic  modulus  .It  tabulated  as  a  function  of  temperature  and  a  linear 
Interpolation  scheme  Is  used  to  calculete  e  value  at  sny  Intermediate 
temperature. 

The  plastic  component  can  be  represented  by  a  series  of  multilinear,  stress' 
plastic  strain  curves.  Lech  curve  corresponds  to  a  different  temperature. 

1. Inoar  interpolation  Is  used  to  calculate  the  strain  for  Into  mediate  values 
of  temperature, 

The  thermal  component  can  bn  written  as 
'll!  ‘  i(t)  1 

wher  u 

a  -  thermal  expansion  coefficient  ■  Integrated  thermal  strain 
It* tween  room  temperature  and  temperature  T  divided  by  the  temperature  T. 

A  table  of  temperature-dependent  thermal  expansion  coefficients  Is  then 
compiled  and  an  Interpolation  scheme  used  for  Intermediate  values. 

Cl 


The  creep  strain  can  be  written  In  a  functional  form  which  has  the  following 
properties : 

1.  Approximates  short  time  creep  behavior  in  primary  and  tertiary 
stages. 

2.  Allows  simple  calculation  of  coefficients  using  analytic  as  opposed 
to  hand  reduction  techniques. 

3.  Matches  the  form  of  the  creep  equations  which  are  either  built  in 
or  can  be  added  to  general  purpose  codes. 

The  most  general  form  of  the  relationship  is  given  in  Equation  3. 


ecr  "  Ao  H-exp^t)]  exp(-A2/T)oA3  +  Bo(tBl)>xp(-B2/T}oB3 

(primary)  ("generalized"  secondary) 


where 

A  ,  A^ ,  A2,  A3,  B0.  Br  B2,  and  B3  are  constants. 


(3) 


It  was  found  that  little  was  gained  by  fitting  the  date  to  both  primary  and 
"generalized"  secondary  portions  of  Equation  3.  Thus,  the  "generalized" 
secondary  relationship  was  adopted.  The  data  from  the  curves  in  Figures  20-24 
was  fit  to  Equation  3  using  a  multiple  linear  regression  procedure  as  outlined 
below  and  described  In  References  7  and  8. 


Given  ecr  *  BQtBl  exp  (-B2/T)oB3 

then 

In  ecr.  =  In  (B0)  +  B^  1  n(t )  -  B2/T  +  B-jln  (o) 

Define  for  each  data  point  (1)  in  each  creep  strain  curve 


B0  *  In  (B0),  =  1  n ( 1 1 )  F^  =  -  1/T,,  F,3  =  ln^) 

then  the  following  system  of  algebraic  equations  can  be  solved  for  the 
unknowns , 


(4) 

(5) 

(6) 
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n  j  n  i  n  «  «  'j 

E  -  5  N  +  B,  IF  +  B-  I  FT  +  n  FT 

cr  o  112  1  1 

1-1  1=1  1*1  1=1 


_  ji  _  n  1  **11  n  |  p  n  •>  ., 

E  ecr  F  =  Bo  1  F  +  B1  1  F  F  +  B  E  F  F  +  B3  E  F  F 

1=1  cr  1  0  1=1  1  1  1=1  1  1  2  1.1  1  1  3  1-1  1  1 


h  _  i?  ^  ^  p  o  p  9 

1  e,r  F7  •  z  Fi  =  B,  I  f!  FT  +  B-  E  Ff  F?  +  B,  E  FT  F'T 

1=1  cr  1  0  1=1  1  1  1=1  1  1  2  1=1  1  1  3  1=1  1  1 


n  i  *3  _  n  <j  n  i  n  o  3  n  ^  3 

1  l  N  =  Brt  I  Ff  +  B,  E  F  Ff  +  Bp  E  Ft  F?  +  B,  E  F?  FT 
1=1  cr  1  0  i=1  1  1  1=i  1  1  2  1=1  1  1  3  1=1  1  1 

(no  sum  on  i) 


In  principal,  a  whole  series  of  equations  such  as  Equation  4  could  be  developed, 
each  corresponding  to  a  particular  temperature  and  stress.  In  order  to  cal¬ 
culate  creep  strain  for  Intermediate  values  of  stress  and  temperature,  one 
can  interpolate  between  curves.  This  procedure  would  most  likely  give 
better  results  in  the  simulation  analysis.  A  single  master  equation  can  only 
give  least  square  type  averaged  values  of  creep  strain.  Unfortunately,  many 
codes,  Including  ANSYS,do  not  allow  for  multiple  creep  equation  entry  and 
thus  the  single  equation  approach  was  taken. 


In  order  to  calculate  the  total  strain  at  any  point  in  time,  one  must  add  all 
of  the  Increments  In  strain  up  to  that  point  In  time.  If  c|ctotai  Is  the 
Increment  change  In  strain  between  time  t  and  time  t  +  At,  then  the  total 
strain  at  time  tQ  Is  given  by  Equation  7. 


'total, 


t 

=  /  de 


total  £total  (Initial) 


Theoretically^ one  could  evaluate  Equation  7  exactly,  given  analytical  forms 
of  the  stress  and  temperature  histories.  In  practice  It  is  much  more  con¬ 
venient  to  simply  evaluate  the  total  strain  by  numerically  stepping  through 
time  In  an  Incremental  fashion.  This  Is  the  procedure  used  in  MARC,  ANSYS 
and  CREEPARHS.  In  particular,  one  calculates  the  stress  and  temperature  at 
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a  specific  point  in  time.  This  allows  direct  calculation  of  the  elastic  and 
thermal  strain.  If  the  particular  combination  of  temperature  and  stress 
Imply  a  loading  state  with  respect  to  the  previous  state,  then  one  can  move 
up  the  plastic  stress-strain  curve  and  calculate  the  plastic  strain. 

If  the  new  temperature  and  stress  state  imply  unloading,  then  the  incremental 
change  In  plastic  strain  Is  zero  and  one  unloads  elastically.  The  calculation 
of  the  creep  strain  requires  the  adoption  of  a  particular  hardening  rule. 

That  Is,  one  must  prescribe  how  the  creep  strain  rate  depends  on  the  creep 
strain,  temperature,  etc.  A  "time-hardening"  rule  was  selected  both  for  the 
sake  of  convenience  and  because  this  rule  1-  thought  effective  when  the 
material  Is  subject  to  temperatures  near  melt  (Reference  5).  In  accordance 
with  the  time  hardening  hypothesis,  the  creep  strain  rate  at  any  time  and 
stress/temperature  level  Is  assumed  not  to  depend  upon  the  current  value  of 
creep  strain.  If  the  stress/temperature  Is  changed  from  o-|/T-|  to  c^/T^  at 
time  t2»  the  creep  rate  Is  determined  at  point  (B)  (Figure  36).  Then  the 
creep  rate  Is  found  from  Equation  8. 
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(3) 


! 


Furctlon  (e  ,  T,  o,  t)  •  0 
cr 

Using  Equation  3  It  can  be  shown  (Reference  9)  that  the  creep  strain  rate 
Is  given  by  Equation  9. 


ecr  -  [A0A1  exp(-Ajt)  exp  (-A2/T)  0*3  ♦  BQB^  tV1  exp  (-B2/T)oB3]  (9) 

The  creep  strain  at  any  point  In  time  t^  Is  given  by  Equation  10. 


4r  ■  “  ccr  <*">  *  <»„  '  Vl> 

Note  that  delayed  creep  recovery  effects  are  not  Included  In  the  model 
because  It  Is  thought  that  the  time  scale  of  the  delayed  effects  Is  long 
compared  to  the  time  considered  In  the  currant  problems. 


(10) 


7.3  C00E-MATERIAL  DATA  INPUT 

The  specific  material  property  data  which  was  Input  Into  MARC,  ANSYS  and 
CREEPARHS  consisted  of  the  following: 

1.  Thermal  stress  versus  temperature 

2.  Piecewise  linear  stress-strain  curves  (Figures  37  and  38) 

3.  Master  creep  strain  curve  (Equation  11) 

Aecr  -  BqBj  t(Bl  ■  1jexp  (-B2/T)o  B3  At  (11) 

where 

B  •  1 .975  x  10*6 
o 

B,  -  0.7056 

°2  ■  ’■J73  «  '°4 
B3  -  2.836 

An  excellent  summary  of  existing  creep  material  property  data  can  be  found 
In  Reference  10. 
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„  7.4  ANALYSIS  VERSUS  SIMULATION  TESTS  -  CREEPARHS  CODE 

The  strain  versus  time  history  as  predicted  by  the  CREEPARHS  code  based  upon 
the  use  of  "zero"  time  +  creep  data  is  plotted  against  the  results  from  the 
simulation  experiments  in  Figures  39,  41,  43  and  45.  The  temperature  and 
loading  versus  time  are  shown  at  the  top  of  the  figure. 

A  0.25  second  time  step  was  used  for  the  creep  analysis  in  Figures  39,  41,  43 
and  45.  In  addition,  the  predicted  curves  are  broken  down  into  thermal, 
mechanical  and  creep  components  so  as  to  illustrate  their  relative  contri¬ 
butions.  Similarly,  the  strain  versus  time  histories  based  upon  10  second 
isochronous  data  is  shown  in  Figures  40,  42,  44  and  46. 

The  simplest  load  temperature  history  leads  to  monotonical ly  increasing  strain 
with  time  as  Is  seen  in  Figure  39.  The  analysis  predicts  a  response  which  is 
in  excellent  agreement  with  the  experiment  until  21  seconds  when  tertiary 
creep  leads  to  ultimate  specimen  failure.  A  more  refined  creep  model  is 
obviously  needed  for  this  analysis.  Examination  of  Figure  40  where  10  second 
isochronous  data  is  used  shows  fair  correlation  through  15  seconds.  This  is 
a  result  of  the  fact  that  thermal  strain  is  dominant  until  20  seconds  when 
creep  failure  occurs. 

As  the  load-temperature  history  becomes  more  complicated,  the  demands  on  the 
material  models  become  more  stringent.  The  results  from  such  an  experiment 
are  shown  in  Figure  41.  Again,  the  correlation  is  encouraging.  In  this 
case,  the  creep  is  insignificant  and  a  failure  due  to  essentially  time 
independent  strain  occurs.  Figure  42  shows  the  predicted  response  when  the 
isochronous  data  are  used.  The  analysis  breaks  down  at  11  seconds  because  of 
the  inability  of  the  stress-strain  curves  to  model  the  plastic  strain  at  the 
given  stress  and  temperature. 

A  more  realistic  simulation  of  an  actual  rapid  heating  problem  would  include 
stress  unloading  due  to  material  property  thermal  degradation.  If  one 
assumes  that  the  material  unloads  elastically,  then  any  residual  strain  would 
be  due  to  permanent  plasticity  and  creep.  The  results  presented  in  Figure  43 
,)  illustrate  the  Importance  of  accounting  for  creep  strain  in  a  problem  where 
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TIME  (SCC) 

CORRELATION  -  EXPERIMENT  VS  "CREEPARHS'COOC  -  TEST  73 

"ZERO"  TIME  DATA 


unloading  occur*.  The  mult*  from  analysis  using  the  Isochronous  data  ara 
4  shown  In  Figure  44.  in  this  cata,  It  1*  elair  that  tha  Isochronous  d#*«  ara 

Inappropriate.  Tha  uta  of  such  data  loads  to  a  prediction  of  fallura  at  a 
vary  early  time  whan  In  fact  no  failure  occur*. 

Tha  most  complicated  load-temperature  hlttory  lead*  to  tha  result*  shown  In 
Figure  48,  Obviously,  In  this  exempt*,  the  qualitative  at  wall  si  quantitative 
behavior  of  tha  tensile  specimen  Is  being  accurately  modeled  by  tha  CRMPAAHS 
code.  A  comparison  with  the  results  shown  In  Figure  46  shows  good  correlation 
through  10  seconds.  At  this  time,  the  Isochronous  based  analysis  again 
predicts  failure. 

7.6  ANALYSIS  •  MARC,  ANSYS,  CACIPARHS  COMPARISON 

This  section  will  Include  a  discussion  of  the  models  that  were  used  In  the 
MARC  and  ANSYS  enslysl s  fol lowed  by  a  comparison  of  the  results  predicted  by 
MARC,  ANSYS  end  CRCCPARM5  for  two  of  the  problems  discussed  In  Section  7,4, 

The  MARC  enelysli  was  performed  on  an  IBM  370/16N  computer  and  the  ANSYS 
analysis  war  performed  on  a  COC  Cyber  74  (6600)  computer,  The  very  elmple 
structure!  model  Is  shown  In  Figure  47.  A  plane  stress,  4-node  Isoperimetrk 
element  (No,  3)  was  used  In  MARC  and  a  similar  ?=0  isoparemetrn  solid  element 
(No.  47  was  used  In  ANSYS,  The  elements  were  issumed  clroply  supported  at 
node  1  and  roller  supported  at  node  'l.  The  element  was  loaded  by  a  uniform 
traction  o  end  the  Poisson  ratio  was  set  equal  to  zero. 
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Figure  4/.  ANSYS /MARL  Structure!  Model 
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Material  properties  modeled  as  time-independent  stress-strain  curves,  a  master 
creep  function  and  thermal  coefficient  of  expansion  versus  temperature  as 

previously  described  were  Input  directly  into  ANSYS.  The  MARC  Input  required 
the  development  of  a  user  written  subroutine  for  both  the  creep  and  the 
temperature  dependent  stress-strain  data.  Care  was  taken  to  make  sure  that 
the  data  input  was  equivalent  in  both  codes. 

A  key  element  in  the  analysis  was  formulating  a  strategy  for  determining  the 
time/load/temperature  steps  required  for  accurate  results.  A  0.25  second 
time  step  was  selected  for  the  CREEPARHS  analysis  and  this  was  the  minimum 
time  step  used  in  the  MARC  and  ANSYS  analysis.  The  MARC  and  ANSYS  analyses 
were  performed  by  marching  through  a  series  of  mechanical  and/or  thermal 
load  steps,  each  followed  by  a  creep  step  where  the  temperature  and  stress 
corresponding  to  the  previous  step  Is  held  constant.  A  rough  determination 
of  the  number  and  timing  of  the  steps  was  made  based  upon  a  review  of  the  behavior 
as  seen  In  the  experimental  data.  The  steps  were  fine  tuned  after  reviewing 
the  analysis  results. 

The  AUTOCREEP  option  was  used  In  MARC.  This  Is  a  routine  which  varies  the 
creep  time  Increment  based  upon  tolerances  for  calculated  ratios  of  Incremental 
creep  strain  to  total  strain  and  incremental  change  In  stress  to  total  stress. 

A  similar  procedure  which  Is  built  into  ANSYS  was  also  used. 

A  comparison  of  the  results  as  predicted  by  the  three  codes  Is  shown  In  Figures 
48  and  49.  These  two  examples  were  chosen  because  all  of  the  various  load 
and  temperature  effects  that  are  demonstrated  In  the  other  tests  are  Inherent 
In  these  problems.  It  Is  seen,  that  the  results  are  almost  Identical.  The 
small  variations  are  due  to  different  methods  of  inputlng  discrete  load  and 
creep  steps  and  slight  differences  In  the  input  plasticity  data.  Any  of  the 
three  codes  could  be  used  to  analyze  these  types  of  problems  with  equal 
success.  The  actual  code  selection,  however,  should  be  decided  based  upon 
problem  complexity  and  cost  considerations . 

Table  7  shows  a  comparison  c>e  tne  number  f  load  steps,  the  total  running 
time  and  the  cost  per  run  for  each  of  the  ihroe  codes.  The  difference  in  the 


Figure  49.  Analysis  Comparison  -  MARC/ANSVS/CREFPAKHS  •  TEST  83 


costs  in  ANSYS  and  MARC  Is  mainly  due  to  higher  prescribed  creep  tolerances 
In  ANSYS  and  the  need  In  MARC  to  repeatedly  access  the  user  written  sub¬ 
routines. 


Table  7 

Code  Execution  Comparison 


The  cost  of  performing  this  simple  analysis  suggests  that  one  would  want  to 
use  a  simple  code  when  possible.  An  efficient  solution  strategy  should  be 
developed  before  a  combined  thermal-elastic-plastic-creep  analysis  of  a  com¬ 
plex  structure  is  undertaken. 
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SIMULATION  EXPERIMENTS 
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x,  / cn  cu  chow  the  results  from  the  many  simulation 

The  following  figures  (50  -  61)  show  the  r  Qf  each 

test/CREEPARHS  analysis  that  were  performed.  Inclu  ed 

figure  is  a  sketch  of  the  stress  end  temperature  histories. 
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Figure  58.  Experiment  vs.  "CFtCEPARHS"  code  -  Test  80  "Zero"  Time  Data 


The  following  consists  of  a  definition  of  the  input  used  in  the  CREEPARHS 
code  followed  by  a  listing  of  the  code.  The  code  was  run  interactively  on 
a  CYBER  CDC  computer.  [Input  temperatures  are  in  (°F).] 


Time  Time 

TTF  =  Run  Time 

TINC  =  Time  Step 

TINT  =  Initial  Temperature  (°F) 

THIN  =  Initial  Thermal  Strain 

NS  =  Number  of  temperature  dependent  sets  of  material  properties. 

TEMP(I)  =  Ith  temperature  for  material  property  input. 

ALPH(I)  =  Ith  Coefficient,  of  thermal  expansion 
Elmod(I)  =  Ith  Elastic  modulus 

NPBT  =  Number  of  breakpoints  for  stress-strain  curves 

EPBK(I,J)  =  Jth  strain  breakpoint  for  Ith  temperature. 

STBKP( 1 ,0)  = Jth  stress  breakpoint  for  Ith  temperature 
EYBK(I)  =  Elastic  modulus  for  Ith  temp  =  Elmod  (I) 

FOR(I)  Ith  creep  coefficient  ec  = 

I  «  1  4  A,  I  =  3  ^  A3 

I  =  2  =?  A2  I  =  4  ^  A4 
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r  ?  i.i',  ■  ■ 

:  C  i 


w 


a  ’Kx 

i  E  "  i  ■  1  >J  '  Ef  :  ■  1  l.l  •  1  n  ■  .  j-  1  F  ’  4  •  •  HI  1  >_1  1  .  I'E  '  :  ■  1  *.’■  1  •  El  1 »  3^  v 

C  Of-Vj'l-n  :•  1  ■  j  I :  I  TLn  IM.  :  Id  •  III  •  r  ;.I.|  •  T  :.F  .  TT  1  ■  IFF 

MM-  0  .  -j  \\ 

T  M 1 1  *  -  0 . 0  VS 

I  EC  - 0 

I  C  T  =  0  •  S 

E  P  :•  T  t  i  =  0 . 0 
THSTN*0. 0 
C RiTh=0 . 0 

READ'  5>  1  »  NS»  NBPT  , T I  NT  «  T 1  Ml »  TH I N 

1  FORMAT  •  .31 3* 2F 10. 2, F 15. 7) 

FEAD'5'2;  SIG1N, TEM IN, TTF . T INC  *  TS I »  TSU, T8F *  TT  I  r 

WP  I  TE (6 »  3  >  SIGIM. TCI ,  TSU. TSF , TEM IN. TT I ,  TTF , T I HC . T I M T . T I M 1  •  T *  1 U 

2  FORMAT  1 2F 10. 5-'2F  10.  5  3F10.5* F10.5) 

3  FORMAT  ( 3X*  ♦MAX  STPE  »  F  1*3.  3i  •'  *  3X«  ♦MR/!  STRESS  OCCUR:-  AT 

1  T IME  =  *  *  F 10. 3, v, 3X»  ^STRESS  REDUCED  AT  T  IME  =  * *  F 1 0 . 3*  *  3! . «  •  STF £  3 S 
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5  ;X»  *  INITIAL  THERMAL  STPA IN=*, E 15. 5) 

Do  4  1=1, NS 

READ1 5»  5;  TEMP'  I ) ,  ALPH ■'  I  >.ELMODC  I ) 

WR I  TEC S.7)  I. TEMP (I >  ,  ALFHs I ) »  ELMODC I > 

4  CONTI  HUE  ' 

READ1 5* 17)  <EPBK<1. JS>, JS*1*NBPT> 

DO  44  fj*l»NS 

READ C 5*  IE)  (STBKPCI  J.  in,  IK=1,NBPT> 

EYBK<IJ)*STBKP(lJ»2)/fEPBKa»2) 

14  CONTINUE  -  -  • 

DO  45  IR-l.NS  ’  " 

DO  45  I  T  =  1 »  KEPT 

EPBK <  IR» IT  >=EPBK( 1 « IT) 

WRI TE(6, 12>  I T, EPBK < IP, IT; » STBKPCIR, IT) 

15  CONTINUE 

5  FORMAT ‘Fl5.5,2E15.5’*  n 

16  FORMAT ( SF8.  1 )  ■>, 

l?  F0PMATC8F8.5; 
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1  E15.3»/,3X,*ELASTIC  M0D=*,El5.3> 

i 2  FORMAT <  3X» *BKPT=*, 13’ 2X, *STRAIN=*,E15. 3, 2X, *3T~tSS  =  M El 5. 3  * 

DO  3  K"  1,4  .  ■  - 

READ*  5,  3  >  FORGO  .  .  ..  ,  , 

8  CONTINUE  .  ’  ■  >, 

9  FORMATCF 15. 5)  ’’  '•  V  ' 

WPITE<6*  14)  f  FOP ^ M ,  M=  1 . 4)  '  M 

14  FORMAT*;  3X»  *CPEEP  CONSTANT***  E  15. 5. /,  3X. ♦STRESS  CONST  ANT =* , 

1  FlO. 5«  '»3X,*TEMP  CONSTANT** ,  E  15. 5,  '  ,  3X,  +TIME  CONSTANT**, 

2  FI 0,5) 

l 1  CALL  ST  I N<  T I M , T I NC , SGI N. TM I N . T INI , S I CAV • T I MAV , TEMAV •  I C  T ) 

CALL  TSTNCEYBk  , TEMP, ALPH, STEGP.EPBl  , FOP  -  HE • AL* EL- SCI N. T M I N *  E R S T 1 1 . 
1THSTN.CPSTN, TIM, TIM1. THIN, SI GAV, T I MAV, TEMAV, T INC*  I  EC  - NSGiBFT •  I C  T  ■ 
I  EC  =  IEC+ 1 
CPSTN*0 . 0 

TOST  Uf-£ES-IU+IitSJ  J  ItjCBS-T  M  ■  & 
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2  ♦THERMAL  STRAIN**.  E15. 5,  /,  3M.  ♦CPEEF  STRA  IN***  E 15. 5*  • 

3  3N*  * TOTAL  STRAI N=*»  E 15. 5 > 

IFCTIM  .GE.  TTF)  GO  TO  15 

GO  TO  11  . 

15  CONTINUE 

STOR  '  ,  .  ,  •  '  • 

END 

SUBROUTINE  to  CALCULATE-TIME, STRESS, TEMPEFhTUFE 

SUE  ROUTINE  ST  I  tu  TIM. T INC. SGI H, THIN*  TINT .SIGAV. TIMflV* TEMA'  * ICT 
COMMON  SIGIN. TEMIN.SIG1 .TSI. TSU.TSF.TTI *  TTF 
S I  G 1  =SG  I N 

TEMA=TfilN  ■  v 

T I  MA  =  T  I M  .*•  :  !:i*  -* 

SIGA  =  SGIN  ■“  •  ’•  v  '  .  ■ 

sg  in=0. 0  •• 

T I M=T IM  +  T  INC  ’  '  ' 

IFf.TSI  .GT.  TIM)  SGIN=S!GIN*TIM/TSI  :  ’ “ 

IF(T IM  .GE.  TSI  .AND.  TSU  .GT,  TIM)  SGIN=SIGIN 

IF(T  IM  .GE.  TSU  .AND.  TSF  .GT.  TIM)  GO  TO  20 
GO  T 0  2 1  j  ■  •  /.  i . j  ,  . 

>8  SGIN=SIGIN*<  1.0-<TIM-TSU;/(TSF-fSU>)  VT  v  " 

ICT*  1  ' 

U  IFCTTI  .GT.  TIM/  TM I N=T INT  +  T IM* <  TEMIN-T I  NT )/TT I 
IF< T I M  .GE.  TTI  .AND.  TTF  .GE.  TIM)  TMIN=TEMIH 

SIGA7=CSGIN*SIGA)/2.  0  / “ .  L  .V % »  . 

T I MAS’* CT I M  +TIMA  >/2.  0  »->'-? •  ilV- 

TEMA V*  <  TMI  N+TEMfl  >  /2 . 0  r  '•?  3-*  '  >  - 

RETURN  >*■  •  ■  ‘  ••~r  tr-TiT-H*  iw.  .  .... 

END  •"  --  -  •  '  '■■-/*  -  - 

SUBROUTINE  TO  CALCULATE  STRAIN  :  tSB:- /- -\i : • 

subroutine  t$tn<  eybk .temp, alph.stbkp, ep6r.for.de. AL* EL. SG IN. 

1  TM IN, EPSTN.THSTN, CRSTN, T IM, TIM1 »  THIN, SIGAV, TIMflV, 

2  TEMAV. T INC, I  EC*  NS, NBPT, ICT) 

DIMENSION  EYBK(l),  TEMPO),  ALPH<  1  >  ,  STBKP(  10,  10),  " 

1  EPBK<10, 10) , DE( 10, 10) , AL  C10) , SIGSTR< 10) , FOR< 1 ) » EL  <  1 ) 

COMMON  SIGIN,  TEMIH.SIGl  .  .  .  •  • 

•  NS1-NS-1  •  •••  \ 

IF  (I  EC  .GT.  0)  GO  TO  2 

DO  1  N=1,NS1  ••  .  >z-- '■**-■» o 

AL  <  N  )*(  ALPHC  N+ 1  )-ALPH(N)  )/'<TEMP(N*l  )-TEMP(N>  )  - 

EL<N)*<EYBK<N*1^-EYBK(N))^<TEMP(N+1>-TEMP<N>> 

DO  1  M=1 , NBPT  »  •  .  - 

1  DE<  N,M>  =  <STBKP<N+1,M)-STBKP(N,M))/'(TEMP(N+1  )-TEMP(N)  ) 

2  DO  3  1*1,  NS1  cr  u  -  i^c  **;  . 

TL  =  TEMP(  I )  ■  ■  •  >.  •  ■  ' 

TU  =  TEMP <  I ♦  1  >  ’"i  -  ' 

NE=  I  . 

IFtTL  .LE.  TMI N  .AND.  TM I N  .LE.  TU>  GO  TO  4 

3  CONTINUE 

4  IF1.  ICT  .GT.  0)  GO  TO  22 
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